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Session  1  -  Monday,  9=30  a.m. 
Chairman  -  A.  Simon 

1 .      On  Numerical  Calculations  of  the  Compression 

of  a  Plasma  Cylinder 

A.  C.  KOLB 
Naval  Research  Laboratory 


2 .       On  the  Containment  for  Infinite  Times  of 

Particles  in  a  Mirror  Machine 
CLIFFORD  S.  GARDNER 
New  York  University 

Post  has  shown  that  a  single  particle,  with  velocity  vec- 
tor not  in  a  certain  escape  cone,  is  trapped  in  a  mirror  machine. 
At  least  such  a  particle  is  adiabatically  trapped:   if  e  (Larmor 
radius)  tends  to  zero,  then  T  (containment  time)  tends  to 
infinity.   This  follows  from  Alfven's  discovery  that  the  magnetic 
moment  v j_/B     is  adiabatically  invariant.   Kruskal's  "all  orders" 
generalization  shows  that  T  becomes  infinite  faster  than  any 
negative  powers  of  e.   Numerical  computations  and  experiments 
have  failed  to  indicate  any  upper  limit  on  T  when  e  is  reason- 
ably small  (neglecting  collisional  effects). 

Prof.  J.  Moser  of  New  York  University  has  recently  proved 
an  important  theorem  in  the  mathematical  theory  of  dynamical 
systems  which  is  applicable  to  this  problem.   It  follows  that 
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for  small  enough  Larmor  radius,  substantially  all  the  slngle- 
partlcle  trajectories  which  are  adlabatlcally  trapped  are  In 
fact  rigorously  trapped  for  all  time. 


5.        The  Boundary  Layer  In  the  Confinement  of  a 

One-Dimensional  Plasma 

W.  PASKIEVICI,  A.  SESTERO,  and  H.  WEITZNER 
New  York  University 

In  this  paper  \ie   describe  attempts  to  obtain  a  simple 
configuration  in  which  a  plasma  composed  of  both  ions  and 
electrons  is  contained  by  static  electric  and  magnetic 
fields.   We  consider  the  simplest  non-trivial  situation  in 
which  the  one-component  fields  depend  on  one  space  coordinate 
only  and  the  particle  acceleration  vector  has  only  the  non- 
vanishing  components.   We  assume  the  contained  plasma  to  be 
at  rest  and  in  thermal  equilibrium  vjith  equal  ion  and  electron 
temperatures  and  with  no  fields  present.   For  the  larger  por- 
tion of  the  paper  we  require  that  there  are  no  trapped  particles 
in  the  boundary  layer  between  the  plasma  and  the  vacuum:  that 
is,  vje  demand  that  ail  particle  orbits  start  and  end  in  the 
contained  plasma.   We  assume  that  the  particles  interact  with 
each  other  only  through  the  average  field  produced  by  all  the 
particles,  that  is,  we  employ  the  Vlasov  equation  for  both 
species. 

We  discuss  the  "charge  neutral"  approximation  and  then 
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give  certain  expansions  useful  In  the  numerical  treatment  of 
the  approximation  and  the  exact  solution.   The  essential 
conclusion  of  a  numerical  computation  of  the  charge  neutral 
curve  Is  that  the  boundary  layer  without  trapped  particles 
cannot  be  essentially  charge  neutral  throughout  the  layer. 
Prom  a  numerical  solution  to  the  full  problem  we  find  that 
the  only  solution  Is  physically  unacceptable  as  It  contains 
a  runaway  Ion  jet  and  large  electric  fields.   We  are  thus 
forced  to  consider  more  complicated  containment  configurations 
and  we  present  two  simple  schemes,  one  with  different  electron 
and  Ion  temperatures,  and  the  other  with  trapped  particles. 
We  offer  these  examples  only  to  show  that  some  kind  of  con- 
finement Is  Indeed  possible;  we  do  not  propose  them  as 
physically  Interesting  --  or  attainable  --  states. 


4.  Runaways  In  a  "Lorentz"  Plasma 

MARTIN  D.  KRUSKAL  and  IRA  B.  BERNSTEIN 
Princeton  University 

A  spatially  uniform  gas  of  electrons  colliding  (Fokker- 
Planck)  with  Infinitely  massive  fixed  Ions  but  not  with  each 
other  Is  considered.   The  development  and  fate  of  runaway 
electrons  Is  treated  appropriately  for  the  asymptotic  limiting 
case  of  small  Imposed  electric  field  E.   In  the  region  In 
velocity  space  of  the  main  body  of  the  electron  distribution 
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(for  velocities  of  order  of  thermal  velocity),  the  collisions 
make  the  distribution  spherically  symmetric  on  a  finite  time 
scale;  this  spherical  distribution  then  satisfies  a  much 
slower  diffusion  equation  {b/ht   ~  E  )  with  a  boundary  condi- 
tion such  that  electrons  ebb  away  by  escaping  to  "infinite" 

-1/2 
velocity  (running  away).   On  a  velocity  scale  of  order  E 

times  thermal  velocity  (the  scale  of  Immediate  competition 
between  the  electric  acceleration  and  the  collisions)  these 
runaway  particles  reappear  as  if  from  a  source  at  zero  velo- 
city, of  which  the  strength  (i.e.  the  runaway  rate)  varies 
slower  than  the  natural  time  scale  {^/bt  ~   E^^     )    in  the  region, 
so  that  the  distribution  function  is  quasi-stationary;  it  may 
be  determined  by  finding  the  unique  solution  once  and  for  all 
of  an  equation  of  "two-way  heat  equation"  type  with  definite 
boundary  conditions.   The  far  runaway  region,  of  velocities  so 
large  that  collisions  are  apparently  negligible  compared  to 
acceleration,  is  also  examined;  the  cumulative  effect  of  the 
collisions  turns  out  to  be  unbounded,  however,  the  particles 
eventually  running  away  with  a  very  slowly  (logarithmically) 
spreading  normal  distribution  of  velocities  orthogonal  to  E. 
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Theoretical  Description  of  the  Astron  E-Layer 
WARREN  HECKROTTE  and  V.  KELVIN  NEIL 
Lawrence  Radiation  Laboratory,  Llvermore 

The  E-layer  of  the  Astron  thermonuclear  device  Is  a 
cylindrical  sheet  of  relatlvistlc  electrons.   The  electron 
density  Is  to  be  such  that  the  self  field  of  the  layer  in 
its  central  region  is  greater  than  the  imposed  external 
field,  thus  providing  closed  lines  of  magnetic  field  in 
this  central  region.   The  injection  system  and  the  external 
magnetic  field  are  to  be  controlled  so  that  the  radius  of 
the  layer  is  approximately  uniform  along  its  length  and  the 
canonical  angular  momentum  of  all  the  electrons  is  about  the 
same.   Making  the  assumptions  that  the  E-layer  is  of  zero- 
thickness,  constant  radius   and  finite  length,  the  electrons 
relatlvistlc  with  unique  canonical  angular  momentum,  and  the 
injection  rate  very  slow  and  the  layer  charge  neutralized, 
the  build-up  of  the  layer  in  time,  the  magnetic  field  config- 
uration, and  the  electron  distribution  in  phase  space  have 
been  calculated. 

Prom  these  results  we  can  conclude  for  our  model,  that 
the  layer  is  stable  against  the  bunching  instability.   The  layer 
is  stable  with  respect  to  the  B^  instability  until  the  self- 
magnetic  field  is  2Qffo   greater  than  that  necessary  for  field 
reversal,  and  perhaps  somewhat  beyond  this  point.   The  prob- 
lem of  providing  the  time -dependent  external  field  necessary  to 
Insure  the  desired  layer  properties  is  defined  but  not  solved 

at  this  time. 
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Session  2  -  Monday,  2  p.m. 
Chairman  -  Burton  Fried 

6.     Trapping  and  Prolonged  Confinement  of  an  Energetic 
Deuterium  Plasma  In  a  Static  Cusped  Magnetic  Field 

J.  OSHER 
Los  Alamos  Scientific  Laboratory 

1  2 

The  results  of  recent  experiments   to  test  trapping   of 

an  energetic  plasma  from  a  fast  coaxial  plasma  gun  In  a  static 
cusped  magnetic  field  will  be  described.   The  time  history  of 
neutron  signals,  produced  primarily  by  deuterons  escaping  out 
along  the  cusped  field  lines  and  producing  D-D  reactions  with 
deuterons  on  the  chamber  wails.  Is  taken  as  the  primary  Indica- 
tion of  plasma  confinement.   This  signal  Is  observed  to  Include 
a  substantial  exponential  portion  with  an  average  lifetime 
(t-,/  )  varying  roughly  as  the  cube  of  the  cusp  magnetic  field 
from  10  \isec   at  2  kG  up  to  1  msec  at  7-'^  kG.   The  long  t.,  / 
observed  even  with  background  Dp  pressures  of  up  to  5  x  10   mm 
of  Hg  appear  to  require  "burnout"  of  neutral  gas.   These  and 
other  measurements  are  consistent  with  5  to  10  keV  average 
D  energy  plasma  being  confined  at  an  average  Initial  density 
of  2  X  10  /cm   In  a  volume  of  «^ 10   cm  with  an  electron 
temperature  >  I.5  keV.    The  confinement  times  observed 


"""   J.  E.  Osher,  Phys .  Rev.  Letters  8,  305  (1962 ) 
^  J.  L.  Tuck,  Phys.  Rev.  Letters  3,  313  (1959). 
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for  energetic  deuterons  are  much  longer  than  one  would  calcu- 
late    for  the  large  nonadiabatic  ion  orbits  (Larmor  radii 
typically  5  cm)  expected  in  the  Mark  II  cusped  field  (machine 
radius  12.5  cm).   A  preliminary  model  proposed  to  explain  the 
long  energetic  ion  confinement  times  is  to  consider  the 
electrons  as  the  fundamental  particle  trapped  in  an  adiabatic 
(for  them)  cusped  mirror  field.   These  electrons  form  a  nega- 
tive potential  well  (measured  ^  1  kV  negative)  which  would 
reduce  the  more  nonadiabatic  ion  loss  rate  through  the  cusps. 
The  well  depth,  occurring  primarily  across  a  sheath  ^  2  cm 
thickness,  would  then  adjust  itself  until  the  modified  ion  loss 
rate  equals  the  nonadiabatic  electron  loss  rates  from  processes 
such  as  Coulomb  collisions  (i.e.,  '^^^^  *-  'Tgg)-   Qualitative 
predictions  for  this  confinement  model  will  be  described. 


^  H.  Grad,  Atomic  Energy  Commission  Report  NYO-9^96,  1961 
(unpublished) . 

E.  P.  Gray  and  I.  M.  Miller,  Applied  Physics  Laboratory, 
Johns  Hopkins  University,  Silver  Springs,  Maryland,  Report 
CPR-010,  1961  (unpublished). 


Ten  Millisecond  Containment  of  an  Energetic 
Plasma  in  a  Cusped  Magnetic  Field 

DONALD  C.  HAGERMAN 
Los  Alamos  Scientific  Laboratory 

A  hydromagnetic  gun  pulse  Injects  deuterium  plasma  into 
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1  -3 

a  biconlcal  cusped  magnetic  field     (axial  cusp  distance 
160  cm).   Neutrons  (~10  ),  mostly  produced  by  energetic 
(~10  kev)  Ions  reacting  with  deuterium  on  the  walls,  monitor 
the  decay  of  the  plasma.   After  Initial  fluctuations  the 
neutron  flux  decays  In  a  reproducible  manner.   Between  2  and 
~4o  milliseconds  after  Injection  the  neutron  signal  decays 
linearly;  the  time  to  half  neutron  Intensity  varies  linearly 
with  the  magnetic  field  strength  and  for  2000  gauss  reaches 
28  milliseconds.   The  linear  decay  does  not  fall  completely 
to  zero  but  changes  to  a  more  gradual  decline.   If  the  base 
pressure  of  the  system  Is  raised  to  6  x  10~-^  Torr  the  neutron 

signal  Is  exponential  with  e-fold  times  of  several  milliseconds, 

3 
These  times  are  10  greater  than  charge  exchange  times;  hence 

burn-out  of  the  background  gas  must  have  occurred.   Present 

estimates  of  plasma  density  are  of  order  10   lons/cc. 


"^  J.  L.  Tuck,  Phys.  Rev.  Letters  3,  313  (1959)- 

o 

D.  C.  Hagerman,  Salzburg  Conference  on  Plasma  Physics  (1961) 

J.  E.  Osher,  Phys.  Rev.  Letters  (to  be  published  In  the 

April  15,  1962  Issue). 
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8.     Comments  on  the  New  Results  in  Plasma  Confinement 

J.  L.  TUCK 
Los  Alamos  Scientific  Laboratory 

The  enclosed  preprint  reports  on  the  achievement  of 
remarkably  long  (and  apparently  stable)  confinement  in  a 
cusped  geometry  filled  by  trapping  from  an  externally  inci- 
dent plasma  jet  (entropy  trapping). 

Previously,  capture  and  confinement  in  picket  fences  has 
not  been  striking.   The  change  seems  to  be  the  result  of 
exceeding  some  threshold  in  the  gun  performance,  presumably 
in  1/2  pv  .   In  the  light  of  these  results  with  picket  fence 
Mark  IIB,  given  in  the  preprint,  efforts  were  redoubled  by 
Hagerman  to  get  the  same  effect  in  the  already  existing  larger 
picket  fence,  Mark  III,  which  has  a  diameter  of  I.5  meters  with 
magnetic  field  2000  gauss.   Now  this  also  has  reached  the  new 
regime.   There  are  some  important  differences  in  performance^ 
however.   The  ratio  Larmor  radius/ machine  radius  is  about  I/5 
in  both  cases.   But  Hagerman 's  confinement  is  much  longer  -- 
10  milliseconds  to  half  Intensity  --  and  the  decay  is,  after 
2  milliseconds,  linear.   When  the  neutron  intensity  has  fallen 
to  -0.1  of  its  value  at  2  milliseconds,  it  then  becomes  exponen- 
tial.  The  dependences  of  the  confinement  times   t   in  these 
differing  regimes  on  magnetic  field  B  are  different,  thus 
T  l/e  oc  B^  in  IIB  (200-2000  M-seg).  while  '^^/2   oc  B  in  III 
(2000  M-sec  -  ^0  milliseconds). 
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One  thing  Is  worth  noting  about  the  10  millisecond  result; 
it  seems  beyond  all  question  to  require  the  achievement  of  burn- 
our.   Whereas  for  the  IIB  machine,  one  might  sneak  out  of  It  by 
stretching  the  numbers. 

It  seems  to  me  that  these  systems  should  snap  back  out  of 
the  burnt  out  state  at  some  point  In  the  decay  when  the  neutral 
gas  flux  from  the  walls  exceeds  the  burnout  power  of  the  decaying 
plasmao   In  this  connection  the  time  taken  for  the  walls  to  return 
the  system  to  l/e  of  the  original  pressure  Is  of  Interest.   From 
the  volume  and  pumping  speed.  It  Is  about  5  sec.   At  this  rate, 
the  contribution  of  background  gas  from  the  gun  Is  probably 
greater.   A  snap  back  effect  has  not  yet  been  observed  with 
certainty,  though  the  change  from  linear  to  exponential  Is  sug- 
gestive of  some  change  In  escape  mechanism  with  density.   It 
should  be  realized  that  the  effect  of  a  sudden  Increase  In  charge 
exchange  Is  complicated  --  for  example,  throwing  particles  out 
randomly  Instead  of  escaping  through  the  cusps,  changing  the 
geometry  of  the  neutron  sources  with  respect  to  the  detectors. 

The  P  In  these  two  machines  are  not  known  and  can  only  be 
estimated  roughly  --  10   In  IIB  and  10  ^  In  III.   Hagerman  has 
fairly  clear  Indications  In  III  that  the  plasma  resides  on  the 
curved  lines  on  the  far  side  away  from  the  gun  (from  the  obser- 
vation that  no  neutrons  are  emitted  from  the  end  walls  on  the 
Inlet  side  of  the  mid-plane).   Such  a  geometry  seems  undlstln- 
gulshable  In  principle  from  that  of  a  mirror  machine  albeit  a 
distorted  one.   Now  the  puzzling  thing  Is  that  the  confinement 
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times  we  observe  are  much  longer  than  those  for  mirror 
machines  at  comparable  p.   For,  as  found  originally  by 
Golovln  and  Joffe  and  now  acknowledged  pretty  well  unlver- 
sally,  mirror  machines  at  P  >  10   seem  to  be  flute  unstable. 
Incidentally  the  flute  Instability  usually  has  the  property 
of  ending  confinement  suddenly.   The  gradually  declining 
signal  seen  In  these  picket  fence  experiments  is  new  and 
suggestive  of  stability.   Can  it  be  that  the  peculiar  warped 
shape  of  the  mirror  machine  existing  in  a  low  P  picket  fence 
imparts  flute  stability  by  some  precessional  effect? 

It  has  been  suggested  by  Osher  that  the  excessively  long 
confinement  time  is  due  to  the  presence  of  a  confining  electro- 
static field,  and  indeed  the  interior  of  the  picket  fence  is 
observed  to  be  1  kv  negative.   He  points  out  that  such  a  field 
is  to  be  expected  since  the  electrons  are  highly  adiabatic;  but 
having  such  large  ion  Larmor  orbits  (~l/5  radius  of  the  machine) 
the  ions  are  not.   The  remarks  of  Berkowitz  et  al .  in  the  Geneva 
Conference  reports  are  noteworthy  on  this  subject. 

One  wonders  how  1  kv  small  potential  could  have  any  effect 
on  5-20  kev  ions.   However, it  may  be  that  the  ion  motion  at 
large  radius  in  the  ring  cusp  is  mainly  tangential  so  that  only 
small  energy  is  available  for  escape  in  the  radial  direction. 
A  small  potential  could  then  turn  such  ions  back. 

It  may  be  worth  observing  that  though  for  a  mirror  machine 
a  bundle  of  planar  coaxial  noncolliding  orbits  is  possible 
(exactly  like  particles  coasting  in  a  cyclotron  which  latter 
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is  after  all  a  truncated  mirror  machine),  no  such  simply 
periodic  noncolliding  orbits  seem  possible  in  the  cusped  sys- 
tem.  Hence  the  long  confinements  observed  are  all  the  more 
interesting,  in  that  they  must  involve  complex  three  dimensional 
intersecting  orbital  motions. 

On  the  subject  of  the  relative  contributions  of  wall 
reactions  and  volume  thermonuclear  reactions  to  the  total 

observed  neutro   signal:   The  wall  reaction  yield  Y   is  given 

a 


by 


Y_  =  n  V  F(T)  (where  F  =  adsorbed  target 

a 


yield  per  deuteron,  T  = 

temperature,  V  volume , 

and  the  thermonuclear  yield  by     (t  confinement  time,  aV 

^Maxwell  averaged  cross 
section. 


Y,  =  |-  n^v(av)  t 


and  the  ratio  Y,/Y  oc  nx.   Calculations  suggest  that  Y,/Y 

b'  a  '^'^         b'  a 

should  be  of  order  of  magnitude  l/lO  in  these  experiments, 
though  Y,  has  not  yet  been  detected  experimentally.   But  as 
density  and  time  are  increased,  so  also  must  the  proportion 
which  is  thermonuclear. 

Which  brings  up  the  constantly  recurring  performance  cri- 
terion -  nx.   Thus  the  Lawson  criterion  for  self  maintenance  in 
DT  is  nr  -   10   (T ~ 20  kev ) =  In  terms  of  nr ,    the  new  performance 
in  picket  fence  is  not  outstanding,  being  comparable  to  those 
achieved  in  the  DCX  and  Table  Top  machines  but  far  below  that 
in  the  Scylla  family.   Indication  of  an  even  higher  nx  in  a 
small  fraction  of  a  mm  in  the  dynamic  pinch  of  Andrianov, 
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N.  V.  Phlllipov^et  al . ,  were  reported  at  Salzburg.   Here  the 
confinement  Is  essentially  Inertlal,  and  we  are  in  a  different 
regime  of  reactors  (I  have  always  called  it  the  Columtus  X 
regime).  Intermediate  between  bombs  and  bottles. 


n(cm   ) 

T(sec ) 

nT 

T.  or  E  kev 

,  _>,    ■  -  ^ . ,  —  -    I 

Neutrons 

DCX 

10^ 

37 

4x10^ 

-300  kev 

--- 

Toy  Top 

lo" 

1.4  xlO 

1.4x  10^ 

~3»5  kev 

~lo5 

Picket  Fence 

I  IB 

10^2 

10-^ 

109 

-5-10  kev 

-10^ 

III 

lo" 

10-2 

109 

■-5 -10  kev 

,.10^ 

Scylla  III 

lo" 

2x10"^ 

2  X  10^^ 

~lo9  kev 

5.10^ 

Andrlanov,  N.V. 

Phllllpov,  et  al 
j 

lo" 

2  X  10"^ 

12 
2x  10^ 

-1  kev 

5x  10^ 

Another  way  of  rating  confinement  performance,  which  amounts 
to  the  same  thing  but  which  is  physically  more  illuminating^ is 
to  express  it  in  terms  of  ion-ion  collision  tim.es,  temperature 
being  assumed  adequate.   For  a  thermonuclear  reactor,  confine- 
ment must  be  for  100  to  1000  ion-ion  collision  tim.es.   (Non- 
thermonuclear  devices  may  not  need  so  many. ) 

In  my  opinion,  such  criteria  as  these  have  been  unduly  neg- 
lected in  setting  out  achievements  in  plasma  confinement »   The 
adiabatlc  compression  family  of  devices  (Scylla,  Pharos,  the  GE, 
the  UK,  the  Aachen  devices,  and  Carlde)  have  in  some  cases  held 
a  plasma  together  for  2  90°  ion  collisions  per  ion  --  at  a 
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temperature  of  -2  kev.   This  Is  not  quite  hot  enough  for  a 
reactor,  or  anywhere  near  long  enough,  but  moving  In  the  right 
direction.   Expressed  In  these  terms,  the  confinements  In  the 
picket  fence  and  any  other  bottles  are  an  order  of  magnitude 
less.   Since  a  hydromagnetlc  Instability  can  end  confinement 
In  less  than  a  collision  time,  observations  of  long  time 
colllslonless  confinement  are  encouraging  to  the  extent  that 
the  confinement  Is  better  than  It  might  have  been. 

Which  brings  us  to  the  concluding  point.   The  plasma 
temperature  in  these  picket  fences  is  probably  already  suffi- 
cient for  reactor  purposes,  though  barely  so.   It  is  probably 
over-optimistic  to  expect  that  these  newly  observed  confine- 
ment times  in  picket  fences  and  their  B  dependence  will  continue 
on  to  higher  densities  and  fields.   But  if  it  should,  then  the 
energy  breakeven  point  for  a  picket  fence  thermonuclear 
reactor  turns  out  to  be  in  the  region  of  n  ~  10  ^,  t  ~  l/4  sec, 
T  ~  10  kev,  and  B  ~  50  kilogauss. 


9«  On  The  Theoretical  Interpretation  of 

Cusped  Containment  Experiments 
HAROLD  C-RAD 
New  York  University 

The  estimation  of  containment  in  a  Picket  Pence  or  Cusped 
Geometry  is  complicated  by  many  factors.   There  are  adiabatic 
and  nonadiabatic  parts  of  both  the  ion  and  electron  populatic^n, 
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each  with  its  own  containment  mechanism  and  nominal  contain- 
ment time.   These  times  differ  by  many  orders  of  magnitude 
and  must  be  balanced  against  one  another  to  obtain  the  actual 
containment  time.    These  Individual  loss  rates  are  strongly 
affected  by  electric  fields  perpendicular  to  B  on  which  depends 
the  type  of  adlabatlc  Invariant  which  governs  reflection  at  a 
mirror  (in  other  words,  the  effective  "hole"  size).    To  be 
effective,  this  potential  does  not  have  to  be  comparable  to  kT 
In  magnitude.   A  different  effect  Is  the  equalization  of  nominal 
electron  and  Ion  losses  which  Is  effected  by  an  electric  field 
along  B;  as  a  rough  rule,  the  net  loss  rate  Is  the  smaller  of 
the  two.   Scaling  of  the  adlabatlc  and  nonadlabatlc  component 
losses  (which  combine  to  give  the  total  nominal  electron  or  Ion 

rate)  Is  radically  different  with  respect  to  n,  T,  B,  as  well 

3 
as  shape  and  size. 

Applying  this  theory  to  the  Osher  and  Hagerman  experiments 

Indicates  that  electrons  are  almost  entirely  adlabatlc  and  Ions 

mostly  nonadlabatlc.   The  observed  thin  sheath  (perpendicular 

component  of  E)  and  corresponding  thin  "hole"  size  (which  was 

predicted  as  a  necessary  requirement  for  good  containment  in 


For  a  general  survey  see  H.  Grad,  "Containment  in  Cusped 
Plasma  Systems",  NYO-9^96,  March  1961. 

^  H.  Grad,  "General  Adlabatlc  Theory",  TID-7520,  p.  402, 
Sept.  1956. 

^  See  H.  Grad,  "Developments  In  the  Theory  of  Cusped  Geometries", 
TID-7520,  p.l48,  Sept.  1956;   more  recently,  "Thermonuclear 
Plasma  Containment  In  Open-Ended  Systems,  NyO-9355^  Sept. i960. 
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4 
the  thermonuclear  range  )  indicates  a  nonadiabatic  ion  con- 
tainment time  which,  although  5~10  times  higher  than  the 
uncorrected  value  (ignoring  the  sheath  potential),  is  still 
much  shorter  than  the  nominal  electron  time »   On  the  other 
hand,  any  adiabatic  ions  will  have  a  containment  time  which 
is  much  longer  than  the  nominal  electron  timeo   Therefore  the 
net  ion  time  may  be  extremely  sensitive  to  the  presence  of  a 
small  number  of  adiabatic  ions,  and  the  range  of  possibilities 
surrounds  the  nominal  electron  time.   Since  the  shorter  of  the 
two  times,  net  electron  and  ion,  is  the  controlling  one,  the 
observed  containment  will  be  sensitive  to  a  change  in  the  net 
ion  time  only  if  this  is  the  shorter  of  the  two;  otherwise 
it  will  depend  almost  entirely  on  electron  containment.   Either 
the  nominal  electron  time  or  a  time  somewhere  between  the 
extremes  of  the  two  kinds  of  ions  is  not  inconsistent  with  the 
observed  decay  times.   However,  the  unusual  scaling  as  B  in 
Osher's  experiment  (which  is  inconsistent  with  either  purely 
adiabatic  or  purely  nonadiabatic  containment)  as  well  as  the 
complicated  transitions  observed  with  time  Indicate  that  the 
experiment  lies  in  a  transition  regime . (transferring  from  non- 
adiabatic to  adiabatic  as  B  increases)  or  else,  as  suggested 
by  Osher,  the  scaling  may  be  due  to  "gun"  effects  which  lead  to 
variations  in  the  state  of  the  injected  electron  plasma. 

In  any  event,  one  cannot  conclude,  at  this  time,  that  there 


IT 

H«  Grad,  "Theory  of  the  Cuspidor",  TID~7505,  Po319,  Feb.  1956. 


-  19 


Is  any  discrepancy  between  theory  and  experiment.   The  most 
likely  method  of  reducing  the  ambiguity  of  the  comparison 
between  theory  and  experiment  would  seem  to  be  to  experim.en- 
tally  determine  the  split  between  adiabatic  and  nonadiabatic 

components,  e.g.,  by  varying  the  relative  field  strengths 

5 
at  the  two  mirrors. 

While  the  observed  containment  time  is  very  interesting 
in  itself,  the  ultimate  significance  is  very  hard  to  predict 
because  of  the  transitional  nature  of  the  present  experiments 
with  regard  to  scaling.   To  become  Interesting  for  thermo- 
nuclear purposes,  one  must  ultimately  reach  a  regim.e  in  which 
the  net  ion  containment  is  much  better  than  adiabatic  electron 
confinement,  and  the  present  experiments  can  only  begin  to 
shed  light  on  this  region. 


-^  See  footnote  1 


10.  R.  F.  Enhancement  of  Single  Particle 

Confinement  in  Cusped  Geometries 
ERNEST  P.  GRAY 
Applied  Physics  Laboratory,  The  Johns  Hopkins  University 

The  possibility  of  plugging  the  ring  and  point  cuspi::,  of 
a  magnetic  monocusp  by  appropriate  r.f.  fields  has  been 
studied.   By  a  numerical  integration  of  single-particle 


This  work  was  supported  by  the  Bureau  of  Naval  Weapons, 
Department  of  the  Navy. 
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trajectories,  the  critical  r»f.  field  strength  for  particle 
reflection  at  the  cusp  --  defined  as  that  field  strength  which 
causes  the  particle  to  be  reflected  just  at  the  cusp  --  has 
been  obtained  as  a  function  of  the  various  parameters,  e.g., 
frequency,  static  magnetic  field  strength,  particle  energy. 
Initial  position  and  Initial  heading.   It  was  found  that 
greatly  enhanced  confinement  can  be  achieved  at  a  relatively 
modest  r.f.  power  level  by  using  r.f.  fields  with  the  electric 
field  direction  more  or  less  perpendicular  to  the  magnetic 
field  lines  of  the  cusp  to  be  plugged,  and  at  such  a  frequency 
that  the  particle  to  be  trapped  will  experience  cyclotron 
resonance  on  Its  way  to  the  cusp.   As  a  typical  example,  a 
proton  of  Initial  energy  1  Kev,  starting  at  the  center  of  a 
monocusp  20  cm  In  radius  and  20  cm  In  length,  with  the  field 
rising  to  30  kllogauss  at  both  ring  and  point  cusps,  cannot 
escape  In  a  single  transit.  If  an  appropriate  30-megacycle 
field  Is  applied  with  a  maximum  electric  field  strength  of 
about  1200  volts/cmo   If  a  Q  of  100  can  be  attained,  no  more 
than  75  kilowatts  of  r.f.  power  are  required,  corresponding  to 
a  power  dissipation  density  of  10  watts/cm   (or  about  10^  of 
the  desired  fusion  power  generation  density).   Although  this 
represents  a  sizable  power  dissipation  density,  it  need  be 
applied  only  over  a  relatively  small  volume  near  the  cusps. 
The  enhanced  confinement  so  obtained  is  not  perfect,  however, 
because  the  particle  eventually  gains  too  much  energy  to  remain 
confined.   After  some  10  to  50  transits  (the  exact  number 
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depending  sensitively  on  the  r.f.  field  strength  and  freq-aency), 
the  particle  escapes  at  an  energy  of  some  20  to  ^0  times  its 
initial  value. 


11.      Some  Problems  of  Nonadiabatic  Plasma  Dynamics 

GEORGE  SCHMIDT 
Stevens  Institute  of  Technology 

Particles  in  adiabatic  motion  are  attached  to  flux  sur- 
faces in  fields  of  axial  symjnetry.   It  is  shown  that  in  the 
case  of  "fast"  changes  in  the  field  (as  measured  in  the  particle 
coordinate  system),  the  particle  tends  to  stick  to  its  flux 
surface,  even  at  the  price  of  sacrificing  the  constancy  of  (x^. 
There  are  circumstances,  however,  when  this  is  not  possible. 
If,  e.g.,  a  quasi -uniform  magnetic  field  is  reversed  in  time 
(second  half  cycle  in  a  9-pinch)  the  first  half  cycle  flux 
surfaces  are  nowhere  to  be  found  in  the  second  half  cycle.   The 
sam.e  is  true  for  a  particle  passing  through  the  median  plane  of 
a  cusp  geometry.   In  these  cases  the  particle  motion  exhibits 
various  peculiarities.   Its  angular  momentum,  e.g., 
mv^n  =  p   -  qrA^  is  strongly  asymmetrical  for  the  change  of  sign 

U  a  a 

of  Aa>    (note  that  Pg  =  const.).   The  lower  bound  of  the  energy 
picked  up  by  a  particle  in  the  second  half  cycle  of  a  0-pinch 
is  calculated  and  found  to  be  very  large  for  electrons  (several 

MeV). 

With  the  breakdown  of  adlabaticity  for  single  particles  the 
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basis  for  hydromagnetlc  treatment  also  breaks  down  for  a 
colllslonless  plasma.   A  simple  self-consistent  field  reversal 
problem  is  investigated.   At  t  =  0  a  cylindrical  shell  of  cold 
particles  is  located  in  a  uniform  magnetic  field  of 
(Pq)  =  n6(P^  -  Pg).   After  the  (slow)  reversal  of  the  external 
field  the  shell  is  found  in  fast  gyration,  setting  up  currents 
weakening  or  even  reversing  the  external  field  insMe.   In 
short,  an  E  layer  is  automatically  established. 

Experiments  performed  and  planned  are  also  shortly 
discussed. 
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12.         Instabilities  In  Inhomogeneous  Plasmas 

E.  A.  FRIEMAN  and  J.  M.  GREENE 
Princeton  University 

We  will  give  a  progress  report  on  our  work  on  the 
stability  of  inhomogeneous  plasmas.   Our  equilibrium  consists 

of  a  magnetic  field  in  the  z-dlrection  supporting  a  plasma 

2 
whose  pressure  is  comparable  to  B  ,  and  whose  temperature 

and  density  vary  in  the  x-directlon.   The  Larmor  orbit  diam- 
eter of  the  particles  is  taken  to  be  much  smaller  than  the 
scale  length  of  the  density  and  temperature  variation.   We 
have  been  considering  wave  propagation  perpendicular  to  the 
magnetic  field.   We  have  shown  that  the  magnetoacoustlc  wave 
is  stable  and  have  identified  an  unstable  mode  with  a  trans- 
verse electric  field.   The  latter  mode  was  previously  reported 
by  Krall  and  Rosenbluth. 


13.  Crossed-Stream  Instability 

E.  G.  HARRIS 
University  of  Tennessee  and  Oak  Ridge  National  Laboratory 

We  consider  a  beam  of  energetic  ions  moving  perpendicular 
to  a  magnetic  field  through  a  background  plasma.  We  show  that 
an  instability  exists  which  Is  similar  to  the  familiar  two- 
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stream  Instability.   If  differs  from  It  In  that  the  electron 
motion  is  essentially  parallel  to  the  magnetic  field  and, 
hence,  perpendicular  to  the  direction  of  ion  streaming.   When 
the  ion  paths  close  on  themselves  and  the  electrons  are  cold, 
we  recover  previous  results.    Our  present  work  includes  the 
effect  of  electron  temperature.   A  nonzero  electron  tempera- 
ture increases  the  instability. 


^      P.  Burt  and  E.G.  Harris,  Phys,  Fluids  ^,  l4l2  (1961). 


l4.      Instability  of  a  Low  Pressure  Plasma  Column 
in  a  Longitudinal  Magnetic  Field 
N.  A.  KRALL  and  A.  SIMON 
General  Atomic 

The  stability  of  a  weakly  ionized,  low  density  plasma  in 
a  uniform  magnetic  field  was  considered  in  a  previous  publi- 
cation.   The  plasma  is  a  finite  conducting  chamber;  there 
is  diffusion  across  the  field  lines  and  streaming  to  the  end 
walls.   An  unstable  helical  mode  results  for  magnetic  fields 
above  some  critical  value. 

The  expression  for  the  critical  field  value  obtained  above 
was  based  on  the  assumption  that  the  eigenfunction  could  be 


* 

This  work  was  done  under  a  joint  General  Atomic -Texas  Atomic 
Energy  Research  Foundation  program  on  controlled  thermo- 
nuclear reactions. 

G.  Guest  and  A.  Simon,  Phys.  Fluids  5,  503  (1962).   See  also 
GA-2671. 
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chosen  to  be  entirely  real.  We  have  now  investigated  the 
validity  of  this  approximation  by  direct  numerical  solution 
of  the  complex  eigenvalue  equations.  Good  agreement  between 
the  numerical  eigenvalues  and  the  analytic  theory  resulted. 
The  variation  of  the  eigenvalues  with  the  parameters  of  the 
problem  was  in  accord  with  the  analytic  theory.  The  nature 
of  the  actual  eigenfunctions  will  be  discussed. 

In  addition  to  the  numerical  check,  one  can  solve  the 
complex  equations  exactly  in  the  limit  of  a  small  parameter 
(which  is  proportional  to  the  magnetic  field).   This  exact 
result  agreed  very  well  in  this  limit  with  both  the  numerical 
results  and  the  approximate  theory  yielding  the  real  and 
imaginary  part  of  the  eigenfrequency  cu  as  well  as  the  eigen- 
functions . 


15.     Dispersion  Relation  for  Fokker-Planck  Equation 

and  the  Effect  of  Collisions  on  Ion  Wave  Instability 
C.  S.  SHEN  and  R.  M.  KULSRUD 
Princeton  University 

The  dispersion  relation  for  electrostatic  oscillations  in 
a  small  external  electric  field  is  derived  on  the  basis  of 
the  Fokker-Planck  equation.   The  effect  of  the  external  field 
and  collisions  on  the  growth  rate  of  ion  wave  instability  is 
discussed  under  the  following  assumptions:   (l)  the  influence 
of  external  field  and  collisions  on  the  perturbed  distribution 
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Is  small;  (2)  the  undisturbed  distribution  functions  of  ions 
and  electrons  are  Maxwelllan  and  shifted  Maxwelllan  respec- 
tively; (3)  the  thermal  energy  of  ions  is  small  compared  to 
that  of  electrons.   The  critical  current  for  the  onset  of 
instability  is  obtained  by  taking  collisions  and  external 
fields  into  consideration.   It  is  found  that  the  external 
field  tends  to  increase  the  growth  rate  of  instability,  but 
the  collisions  tend  to  damp  it. 


16.  Secondary  Plasma  Growth 

T.  K.  FOWLER 
Oak  Ridge  National  Laboratory 

It  is  pointed  out  that,  in  a  plasma  containing  electrons 
sufficiently  energetic  to  ionize  neutrals,  the  cold  ion  density 
increases  abruptly  above  a  critical  neutral  pressure  such  that 
the  ionization  time  is  le  ss  than  the  time  for  cold  ions  to 
drift  out  of  the  plasma.   The  ion  density  grows  until  the 
accompanying  hot  electrons  reduce  the  ionization  time  by  clean- 
ing out  the  neutral  population,  a  "secondary"  burnout.   This 
phenomenon  appears  to  occur  in  DCX-1  and  possibly  accounts  for 
a  puzzling  discontinuous  change  in  density  of  a  plasma  formed 
by  electron  cyclotron  heating  when  the  pressure  exceeds  a  criti- 
cal value.   In  the  latter  case,  the  plasma  potential  would  be 
predicted  to  change  sign  when  the  density  jump  occurs. 

-  27  - 


Session  4  -  Tuesday,  2  p.nio 
Chairman  ~  Edward  Harris 

17-     Q-uantum  Mechanical  Kinetic  Equation  of  a  Plasma 

BURTON  D.  FRIED  and  H.  W,  WYLD,  Jr, 

Thompson  Ramo  Wooldrldge  Inc. 
Canoga  Park,  California 

The  method  of  deriving  kinetic  equations  for  a  classical 
plasma  developed  by  Rosenbluth,  Rostoker,  Lenard  e_t  aJ.  ,  is 
generalized  to  the  case  of  a  plasma  so  dense  that  the  effects 
of  the  uncertainty  principle  and  the  Paull  principle  become 
significant.   In  place  of  the  BGB'ZK  hierarchy  of  equations  for 
the  distribution  functions,  f  ,  one  has  a  similar  coupled  set 
of  equations  for  Wigner  distribution  functions,  or  (more 
convienently  it  turns  out)  for  the  Poijrier  transforms  of  these. 
Since,  in  general,  there  is  no  small  expansion  param.eter  analo- 
gous  to  the   nAr^   of  the  classical  case,  neglecting  the  three 
particle  correlations  has  not  the  same  justification  as  in  the 
classlc>al  case.   Nonetheless,  doing  so  leads  to  a  pair  of 
coupled  equations  from  which  one  can  deduce  m.ost  of  the  informa- 
tion and  results  of  other  fashionable  approximation  methods. 
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l8.        Kinetic  Theory  of  an  Electron-Photon  Gas 

He  DREICER 
Los  Alamos  Scientific  Laboratory 

A  Bolt zmann -like  transport  equation  for  the  distribution 
of  electron  moments  is  derived  for  an  electron-photon  gas  in 
which  electrons  emit,  and  absorb  photons  Independently.   Under 
equilibrium  conditions  the  equation  is  satisfied  by  a  Maxwellian 
distribution  for  the  particles  provided  the  frequency  distribu- 
tion of  the  photons  follows  the  Planck  law  at  the  same  tempera- 
ture.  Application  is  made  to  the  case  of  negative  absorption 
of  cyclotron  radiation  by  the  electrons  of  a  colllsionless 
plasma  situated  in  a  magnetic  field.   Negative  absorption  (or 
raaser  action)  can  occur  when  the  electron  velocity  distribution 
possesses  an  appropriate  positive  slope,  and  provided  transit 
time  broadening  of  the  cyclotron  radiation  is  sufficiently 
small.   Under  these  conditions  the  electrons  radiate  more 
energy  by  induced  emission  than  they  absorb  by  true  absorption, 
and  in  the  process  their  velocity  distribution  relaxes  to  a 
state  closer  to  the  equilibrium  Maxwellian  distribution.   This 
relaxation  process  has  been  followed  with  the  help  of  the 
Boltzmann-like  transport  equation  described  above  by  assuming 
either  a  Gaussian  or  a  Lorentz  frequency  distribution  for  the 
amplified  radiation  field.   The  result  indicates  that  this 
relaxation  process  gives  rise  to  a  broadening  of  the  amplified 
line.   The  rate  at  which  the  velocity  distribution  relaxes  in 
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the  neighborhood  of  its  maximum  positive  slope  Is  proportional 
to  the  Intensity  of  the  radiated  line,  and  Inversely  proportional 
to  the  cube  of  the  line  width  of  the  radiation,  i.e.,  (Aod)". 


19.     Thermalizatlon  of  Fast  Ion  Groups  In  a  Plasma 

HERBERT  C.  KRANZER 
Adelphl  College 

The  process  of  thermalizatlon  of  a  single  fast  Ion  Injected 
Into  a  plasma  in  equilibrium  was  studied  in  (l).   The  procedure 
used  was  numerical  integration  of  the  fully  linearized  Fokker- 
Planck  equation  with  time -independent  coefficients  corresponding 
to  an  undisturbed  plasma.   In  some  proposed  thermonuclear 
reactors,  however,  a  whole  group  of  ions  would  continuously 
transfer  energy  to  the  plasma  during  therm.alizatlon,  signifi- 
cantly altering  the  coefficients  in  the  Pokker-Plance  equation. 
Accordingly,  the  computations  of  (l)  have  been  repeated  with 
two  modifications.   First,  the  initial  ion  distribution  is 
taken  as  Maxwell Ian  rather  than  as  a  delta  function.   Second, 
the  plasma  ions  and  electrons  are  allowed  (separately)  to 
Increase  in  temperature  as  the  injected  ions  cool  down;  the 


Research  supported  by  the  Atomic  Energy  Commission  a.Y\r>   the 
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Kranzer,  H.  C,  "Thermalizatlon  of  a  Past  Ion  in  a  Plasma", 
Physics  of  Fluids  4,  21^-220  (196I). 
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total  energy  of  the  system  remains  constant.   In  general, 
the  mean  time  of  thermallzatlon  for  a  group  of  fast  Ions  Is 
found  to  be  longer  than  for  a  single  Ion  at  the  same  tempera- 
ture . 


20.  Thermal  Relaxation  In  a  One-Species, 

One -Dimensional  Plasma 
JOHN  DAWSON 
Princeton  University 

1  2 

The  theories  of  Lenard  ,  Rostoker  and  Rosenbluth  ,  and 

Balescu  predict  that  stable  non-Maxwell  distributions  do  not 
relax  to  a  Maxwell  distribution  for  a  one-species,  one- 
dlmenslonal  plasma.   These  theories  consider  only  the  two- 
particle  correlation  function  and  apply  only  In  the  limit  of 
a  large  n, ;   n.  Is  the  number  per  Debye  length.   This  prediction 
has  been  verified  by  numerical  experiments  on  a  one -dimensional 
plasma  model.   Systems  containing  1000  to  3000  particles  were 
started  with  square  velocity  distributions «   The  number  of 

particles  In  a  Debye  length  ranged  from  2.5  to  20.   It  was 

2 
found  that  the  time  for  relaxation  to  Maxwelllan  went  as  n, . 


^     A.  Lenard,  Ann.  Phys.  1_0,  390  (1960). 

o 

N.  Rostoker  and  M.  N.  Rosenbluth,  Phys.  Fluids  3^  1  {i960) 

^   R.  Balescu,  Phys.  Fluids  3,  52  (196O). 
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This  relaxation  is  much  slower  than  the  relaxation  of  a  test 
particle  to  the  background  distribution  which  goes  as  n , • 
A  transient  relaxation  also  showed  up  which  required  a  time 
of  the  order  of  l/oj  ,  and  which  pulled  the  corners  of  the  dis- 
tribution  in  a  distance  AV  =  {2va'^D/n)   '       where  a  is  the  charge 
on  a  sheet  and  D  is  the  Debye  length.   The  distribution  also 
develops  a  tail  which  extended  the  distance  AV  beyond  its 
initial  edge.   It  was  also  found  that  though  the  mean  distribu- 
tion function  remained  square  for  quite  a  long  time,  there  were 
rapid  fluctuations  about  this  distribution  and  the  characteristic 
time  for  these  was  of  the  order  of  l/co  . 

J.  Dawson,  Phys .  Fluids  "b Md    (1962). 


21 .  Irreversible  Processes  in  a 

One -Dimensional  Plasma 

0.  ELDRIDGE  and  M.  FEIX^ 

General  Atomic 

We  investigate  the  test  particle  problem,  and  the  evolution 
of  a  non-Maxwellian  velocity  distribution,  using  a  one- 
dimensional  model  of  a  plasma  which  allows  us  to  calculate  with 
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an  electronic  computer  the  exact  trajectories  of  particles. 
We  obtain  from  the  usual  asymptotic  theory  the  Fokker-Planck 
coefficients  for  a  set  of  test  particles,  I.e.,  particles 
which  are  singled  out  from  an  equilibrium  distribution.   The 
numerical  "experiments"  show  that  the  relaxation  time  for  test 
particles  varies  as  nD,  the  number  of  particles  In  a  Debye 
length,  as  predicted  by  the  theory.   The  magnitude  of  the 
coefficients  also  agrees  with  theory  for  4  <  nD  <  20.   The 
first  order  theory  predicts  that  the  detailed  balance  between 
drag  and  diffusion  is  valid  for  all  stable  distribution,  not 
only  at  equilibrium.   So  any  change  in  the  velocity  distribution 
should  be  at  least  a  second  order  effect.   We  study  the  evolu- 
tion of  two  non-Maxwell Ian  distributions,  using  the  entropy  as  a 
measure  of  the  deviation  from  a  Maxwellian .   The  "experiments" 
show  a  large  entropy  change  during  the  first  half  plasma  period 
for  all  values  of  nD.   The  relaxation  then  goes  at  a  much  slower 
rate . 


-  33  - 


NOV  26  1362 

DATE  DUE 

1 

CAYLORO 

PRINTED  IM  U    S    * 

/ 


NYU 

NYO- 

10h2k 


c.l 


New  York  University.    Courant 
Institute   of  Mathematical 
Sciences.       Mflcrnflt-.n-Flnl  ri    Dvnft- 

NYU  c.l 

NY0-10i;2ii 

New  York  University.  Courant 


AUTHOR  Inst,  of  Math.  Sciences 
Magneto-Fluid  Dynamics  Div. 


TiTLESummaries  of  papers  pre- 
sented  at  the  meeting  on 


trollud  nuclear  fusion 


ct&a<r- 


N.  Y.  U.  Courant  Institute  of 
Mathematical  Sciences 

4  Washington  Place 
New  York  3,  N.  Y. 


